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1. Introduction. 

Mechanical systems with constraints on the velocities, called "nonholonomic con- 
straints", have paramount importance in engineering, in particular in robotics, vehicular 
dynamics and motion generation, so they are actively and deeply studied also from the 
theoretical point of view, see the recent Bloch-Krishnaprasad-Marsden-Murray [BK], the 
references therein, and Kupka-Oliva [KO]. We are not going to review the story of these 
studies, which started long ago with the work of Lagrange, let us just mention Chetaev 
among the most important contributors. 

The mathematical model kept essentially unchanged until about 20 years ago when a 
new dynamics of velocity constrained mechanical systems was introduced by Kozlov [K], 
and was reported in the beautiful "Encyclopaedia of Mathematical Sciences" edited by 
Arnol'd, see [A] Ch. 1, §4. This new mechanics was called "vakonomic" being "variational 
axiomatic kind" . 

The paper Lewis-Murray [LM] studies the two mechanics from the theoretical point 
of view, and deals with a ball on a rotating table analytically, numerically, and experimen- 
tally. The experiment supports the nonholonomic framework, moreover the authors say 
at page 808: "we were not able to produce any vakonomic simulations which resembled 
the experimental observations...". However, they also say at page 809: "Certainly, a more 
careful and exhaustive experimental effort on systems other than a ball on the rotating 
tabic would be valuable in providing data which would allow for a fair comparison on the 
nonholomic and vakonomic methods." 

The previous paper Kozlov [K], part III §5, relates the two different theories to the 
different ways the constraints can be realized, e.g. by large viscosity or additional masses, 
and it says "...vakonomic dynamics, which is an internally consistent model that can be 
applied to the description of the motion of any mechanical systems, is as 'true' as traditional 
nonholonomic mechanics. The issue of the choice of model for each particular case is 
ultimately resolved by experiment." 

The aim of the present paper is to suggest the contrary: perhaps vakonomic mechanics 
is not satisfactory as a model for velocity dependent constraints. This opinion is based 
on the main example used by Kozlov to support his dynamics, namely the skate on an 
inclined plane, whose nonholonomic behaviour is rated "paradoxical" in [A] p. 19, and, at 
page 36, to be compared with the vakonomic motion he studies. 

The paradox seems to come from the fact that for the nonholonomic dynamics "...on 
the average the skate does not slide down the inclined plane...". 

In Section 2 we recall the dynamics of natural systems with nonholonomic constraints, 
that is Lagrange equation with multipliers, and we put it in normal form. It is a conservative 
dynamics. Moreover, we see that it is reversible, namely the set of solutions is invariant 
under time inversion. The example of the nonholonomic skate ends the section. 

Section 3 deals with the same systems but with vakonomic dynamics. Now we have 
many more solutions since new "latent variables" A appear ([A] p 37, and the paper [K] 
III at page 44, loosely relate this fact with the "unobservable quantities... in quantum 
mechanics"). This situation, quite strange for classical mechanics, is not overcome by 
passing to the Hamiltonian framework. Anyway, we prefer to keep our discussion in the 
Lagrangian framework, which is, perhaps, clearer for our purposes. 
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We show that also the vakonomic dynamics of natural systems is reversible, at least 
in a suitable sense, and we see no reasons why it should not be reversible, as nonholonomic 
dynamics is. The section ends with the study of some vakonomic solutions to the skate 
which seem paradoxical to me for all initial values of the "latent variable". 

Kharlamov [Kh] considers also the skate on a plane to criticize the vakonomic mechan- 
ics. His plane is horizontal, so he is able to obtain explicit expressions for the vakonomic 
solutions which correspond to the nonholomic uniform motions along straight lines. The 
author concludes that the vakonomic skate has a "fanciful track". However, the uniform 
motions along straight lines are vakonomic motions. This fact does not "save the situa- 
tion" according to Kharlamov since "there is no basis on which to assign a specific value" 
to the latent variable at the initial time. 

Finally, Section 4 deals with more general nonholonomic systems and revisits the 
nonholonomic skate by adding some dissipation on the rotational degree of freedom only. 
Now, the skate slides eventually down. 



2. Nonholonomic dynamics of natural systems. 

The Lagrange equation of Classical Mechanics is the celebrated 

The Lagrangian function L is assumed of class on some open connected subset of 
M X X M^. 

In particular we are interested in the natural Lagrangian functions: 

L{t, q, q) = lq- A{q)q - U{q) , A{q) = A{qf > , (2.2) 

where A{q)'^ is the transpose matrix of the NxN positive definite matrix A{q), the central 
dot is the usual scalar product, ^q-A{q)q is called the kinetic energy, and U (q) the potential 
energy. 

If one asks the mechanical system to obey the nonholonomic constraint equation 

B{q)q = (2.3) 

where B{q) is annxN, with n < N, full rank matrix at each q, then one has a nonholonomic 
system, which we briefly call "nonholonomic natural system", whose dynamics is ruled by 
the constraint (2.3) together with the following classic Lagrange equation with multipliers 

= (2.4) 

The multipliers are the components of the new unknown function /x(t) with values in MJ^. 
In the sequel, we assume S e as L. For natural systems the dynamics is then ruled by 
equations of the following form 

A{q) {q + T{q)[q,q]) +VU{q) = B{qfi, 

(2.5) 

B{q)q = 
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where V is the gradient operator, and T{q) : IR^ x — )■ R-^, (u^v) i-^ r(g)[w, v], is a 
bihnear symmetric map whose components are called Christoffel symbols in the classical 
books. 

We easily check that it is a conservative dynamics, namely the total energy ^q-A{q)q+ 
U{q) is a first integral of (2.5). 

For brevity, let us omit the functional dependences for a while. The N x N matrix A 
is symmetric and positive definite, in particular then invertible, so we can define the new 
n X n matrix C := B A~^ B^ which is nonsingular, indeed = =^ = ^ ■ C ^ ~ 
^ ■ B A-^ B'^^ = B^^ ■ A-^ B^^ =^ B'^i = since A'^ is positive definite as A; finally 
B'^^ = =^ C = since B has full rank. 

It is easy to see that the following {N + n) x {N + n) matrices 



A [A-^-A-^B'^C-^BA-^ -A'^B'^C'^ 

-B h [ -C-^BA-^ -C-i 



(2.6) 



are inverse of each other. Thus the following system in the unknown functions q{t),iJ,{t), 
which is obtained from (2.5) by differentiating the second equation with respect to t, 



-B J U J + l -B'{q)[q,q] i ' " ^''^^ 



is equivalent to 



q \ _ (A-^- A-^B^C-^BA-^ -A-^B^C-^ \ ( -A{q)T{q) [q, q] - VU{q) \ 
-C-'BA-^ -C-i J [ B'{q)[q,q] J ' 

(2.8) 

Reinstating the functional dependences and defining suitable new functions /, g we get the 
following "normal form" which is quadratic in q 

Q = Hq) + 9{q)[q,q] 

(2.9) 

= h{q) + k(q)[q, q] 

We can get rid of the last equation and the unknown 

Since the second equation in (2.7) is obtained from (2.3) by differentiation, we have 
B{q) q = const along the solutions of the first equation in (2.9). So, the constraint equation 
B{q) q = is satisfied provided it holds at some time to- Therefore, we arrive at the 
conclusion that all nonholonomic motions can he obtained from Cauchy problems of the 
following kind where an additional condition restricts the choice of the initial velocity 

(2.10) 

(^(^o) , q{to)) = {qo, qo) G {{u, v) : B{u) v = 0} 

If q{t) is a solution to (2.10), then the function r{t) = q{—t) is also a nonholonomic 
motion which has initial data {r{—to),r{—to)) — {qo, —qo) at time —to- So the set of non- 
holonomic motions is invariant under time reversal, and we can say that the nonholonomic 
dynamics of natural systems is reversible. 
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Example: the nonholonomic skate. Consider a skate (an homogeneous material 
segment) on an inclined plane with Cartesian coordinates x, y. The x-axis points downward 
while the y-axis is horizontal, x and y will denote the coordinates of the center of the 
skate which has another Lagrangian coordinate: the rotation angle 4> h makes with the 
unit vector of the x-axis. We have a natural Lagrangian function as in formula (2.2) 

L{t, X, y, (t>, X, y,4>)^^ (x^ +y^ + <p^^+x. (2.11) 

Suppose the center of the skate can only have velocities parallel to it, namely we 
consider the constraint equation 

X sincj) — y cos(f) = . (2.12) 

The nonholonomic equations (2.5) in the actual case are 

X — 1 = n sin (f) 

y^-l_i cos (j) , . 

= ^'^■^'^^ 

X sin (j) — y cos = 

We easily get the following Cauchy problem which particularizes (2.10) to our system, to 
the initial time to = 0, and to some special initial conditions (considered in [A] p. 19) 

X = cos^ — sin {x cos + y sin 0) 
y = cos sin + cos {x cos + y sin 0) 

0=0 ^ ■ ' 

{x{0), y{0), 0(0), i(0), y(0), 0(0)) = (0, 0, 0, 0, 0, w) 

with w ^ 0. Notice that the initial data above are acceptable since they satisfy the 
constraint at t = 0: x{0) sin 0(0) — y{0) cos 0(0) = 0. The solution is easily found 

x{t) = sin^ wt , y{t) = (wt - \ sin 2wt] , 0(t) = wt . (2.15) 
2w'^ 2w'^ \ 2 / 

It is a cycloid. As remarked in [A] p. 19, on the average the skate does not slide down the 
inclined plane: < x{t) < l/2w'^. 

We are also interested in the motion starting with the conditions (0, 0, 7r/2, 0, j/o, 0), 
which also satisfy the constraint 

x{t) = 0, y{t) = yot, 0(t) = 7r/2. (2.16) 
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3. Vakonomic dynamics of natural systems. 

As in Section 2, we consider the natural Lagrangian L{t,q,q) in (2.2), and the con- 
straint (2.3) with a full rank n x N matrix B{q). As is well known, the nonholonomic 
system (2.5) is a consequence of d'Alembert's principle which is not variational. If we 
adopt a variational approach by requiring the motion to be a stationary curve of the ac- 
tion functional q{-) i— )> J^^^ L{t,q{t),q(t))dt among all curves having the same end points 
and satisfying the nonholomic constraints, then wc get a vakonomic motion (see [A] pp. 
32, 33, 34 and [KO] for details). Moreover, the motion t h- )■ q{t) is vakonomic if and only if 
there exists a smooth curve t i-> X{t) in R", defined on the same time interval, such that 
the pair {q{-), X{-)) is a solution to the (unconstrained) variational problem associated to 
the Lagrangian function 

Lit, q, q, A, A) = Lit, q,q) - X- Biq) q (3.1) 

(the scalar product being in R"^ as the values of the new unknown A(t)). Namely, the 
vakonomic motions can be obtained by the Lagrange equations 



(3.2) 



d dh 


dL 


dt dq 


dq 


d dL 


dL 


dt~dX~ 


dX 



As in the previous section we speak of natural systems, now vakonomic. Their dynamics 
is ruled by 

A(g)(g + r(g)[g,g]) + VC/(g) - i?(g)^ A - A • i?'(g) g + ^ (A • i?(g)g) = 

(3.3) 

S(g)g = 

where, for each i e {!,..., N}, the i-components of two of the previous expressions are 

a=l i=l \ " / i a=l i=l 

As in Section 2 we can derive the second equation in (3.3) with respect to t and write 
the system obtained from (3.3) in this way as 

for a suitable F. Again the iN + n) x iN + n) matrix is invertible since A(g) = A(g)^ > 
and Biq) has full rank at any point g, and the inverse is as in (2.6). So we can solve (3.5) 
with respect to (g. A) and get an equivalent system in normal form 



G'(g,g,A). (3.6) 
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The function B{q)q is a first integral, i.e. constant along the solutions. If B{q{to))q{to) = 
at the initial time to, then we have a vakonomic motion, moreover all vakonomic motions 
can be obtained in this way. However, now A (to) is arbitrary in the Cauchy problem. It 
seems that we have too many vakonomic motions. 

As in the nonholonomic dynamics of natural systems the total energy is conserved. 
However, now we cannot get rid of A. 

In order to discuss the time reversibility, let us consider the new unknown function 
cr{t), with &{t) — X{t). Then Kozlov function (3.1) becomes quadratic in (g, &) 

L{t, q, q, a, cr) = ^q ■ A{q)q - U{q) - a ■ B{q)q (3.7) 
and the related equations 

^ " ° (3.^ 



d dL 


dL 


dt dq 


dq 


d dL 


dL 


dt da 


d^ 



dt da da 

give directly (3.5) with a instead of A, a instead of A, and no dependence on a 

A(q) -B{q)T\ f ij 



-B(,) ) , + ni,i,i) = (3.9) 



This system can again be put in normal form 

9 



^, G(g,g,a) (3.10) 

with the first integral B{q) q. Whenever this first integral vanishes, we get a vakonomic 
motion and all vakonomic motions can be obtained in this way. Finally, the solutions to 
(3.3) are the pairs (g(t),A(t)) = {q{t),a{t)). 

The equation (3.10) is quadratic in the velocities, thus for any solution we have a cor- 
responding reversed one: if (g(t), A(t)) is a vakonomic motion, then so is (g(— t), — A(— t)). 
Perhaps one can argue that these solutions are not necessarily physically acceptable at the 
same time, so accepting one and rejecting the other. I think that this point should be 
somehow justified: without a good reason why should the vakonomic dynamics of natural 
systems miss the feature of time reversibility? Would not the very lack of time reversibility 
be a sufficient reason to reject the vakonomic model? 

Example: the vakonomic skate. Consider again the skate as in Section 2, now 
vakonomic. The above Kozlov function (3.1) for this mechanical system is 

L{t, X, y, (j), X, y, 0, A, A) = ^ + y^ + + a; - A (i sin - y cos 0) . (3.11) 
The vakonomic equations are 



X — X sin (j) — cos = 1 
y + \ cos (j) — X(p sin = 
+ A (i; cos 4> + y sin 0) = 
X sin (j) — y cos = 



(3.12) 
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We easily get the following system 



X = —(j) sin 4> {x cos (/> + y sin cf)) + X<p cos (p + cos^ (f) 
y = +4) cos (f) {x cos 4> + y sin 0) + A sin + sin cos > 
= —A (x cos + y sin 0) 
A = — sin — (a; cos + y sin 0) 



(3.13) 



which has the first integral 

x sine/) — y coscj) = a- (3-14) 

The solutions to (3.13) with a = are all the vakonomic motions. 

The vakonomic motions to be compared to the nonholonomic motions (2.15) (we use 
the plural since there is a parameter w), according to Kozlov have A(0) = (this choice 
is probably justified by the results it gives), and is studied in [K] III, Section 3. Their 
features are also reported in [A] p. 35, 36. Here, let us just quote from there that the skate 
slides monotonically down for t > 0, and almost all solutions tend to turn sideways: (f){t) 
converges to one of the points tt/2 + mn (m G Z), as t ^ +oo. Let us add that for t < 
it slides monotonically up. 

Now, let us go to the nonholonomic motions (2.16). First of all, let us immediately 
remark that they have no corresponding vakonomic motions the extra variable A notwith- 
standing, unless ijQ = 0, the "vakonomic equilibrium", for which the "latent variable" A 
has the precise meaning of reversed time: 

x(t) = 0, y(t) = 0, 0(t)=7r/2, X(t) = X(0)-t. (3.15) 

Consider now (3.13) with the initial conditions 

(x(0), 1/(0), 0(0), i(0), m, m, A(0)) = (0, 0, 7r/2, 0, yo, 0, Aq) , yo . (3.16) 

For the moment assume that Aq 7^ too. Then (3.13) give :r(0) = 0, the second and the 
third derivative of the functions 0(t) and x{t) respectively, at t = 0, are 

0(0) = -Ao yo sin(7r/2) = -Xq yo , x^^\0) = -0(O)|/(O) sin2(7r/2) = Aq (3.17) 

Therefore 

'x{t) = Xoyit^/3\ + o{t^) 
AoT^O =^ { y{t) = yot + o{t) (3.18) 

0(t) =7r/2- AoyotV2 + o(t2) 

If we accept the reversibility of vakonomic dynamics we arrive at the conclusion that either 
the solution we are analyzing, or the one obtained by time reversal, slides monotonically 
up in a full neighbourhood oft = 0. I guess we can rate this behaviour paradoxical given 
the initial conditions above. 
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Next, consider the other possibihty Aq = 0. In this case (3.13) gives 



0(0) = 0, (P^^\0) = -A(0) m sin(7r/2) = yo sin^{n/2) 



x^^^ (0) = 0, a;(^^ (0) = 0, x^^^ (0) = -0^^^ (0) y{0) sin^ (tt /2) = -yl 



(3.19) 



So 



Ao = 



' x{t) = -yo' /4! + o{t^) 
< y{t) =yot + o{t) 



(3.20) 



I (/){t) = 7r/2 + yot^ /3! + o(t3) 

which seems another paradox: the skate slides up for t > (smaU enough) while both x{0) 
and 4){0) vanish, and A(0) = too. 

Of course, also the "body-plus-fluid" , considered in [K2] p. 598 to answer the criticism 
in [Kh], has the same strange behaviour as well as any other physical system we may think 
should obey the same dynamics. 



4. General nonholonomic dynamics. 

For the final example, equations only slightly more general than those considered in 
Section 2 are needed. But full generality can be obtained with only little work from what 
described in Section 2, so we proceed to it. 

The Lagrange equation whith multipliers (2.4) is also studied for general Lagrangian 
functions L{t, g, g), and has been also extended to constraints possibly non-linear in g by 
the work of Chetaev and others. Notice that justifying the extension was a difficult job, see 
the recent Cardin-Favretti [CF] , and Massa-Pagani [MP] ; but concrete non- linear examples 
can be given, see Benenti [B]. 

The general equations for nonholonomic dynamics are 

a^^^'^'^)- 9^^^'^'^)= U^'''''V ^ (4.1) 
b{t, g,g) = 

where, for each t, the values g(t),g(t) G M^, ii{t),b{t,q{t),q{t)) G W^. The standard 
conditions to put (4.1) in normal form are 



det S{t, g, g) ^ where S{t, g, g) := -^(^, ^) at aU (t, g, g) , (4.2) 



and det R{t, g, g) 7^ at all (t, g, g) , 

1 rr. db (4.3) 

where R := D S'^ , with D{t,q,q) := —{t,q,q) ■ 
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Indeed, differentiating the second equation (4.1) with respect to t, and introducing 
suitable new functions r and s, the system (4.1) gives 

S{t,q,q)q-r{t,q,q) = D{t,q,q)'^ /i 

(4.4) 

-D{t,q,q)q- s(t,q,q) = 

Next, we check at once that the following two {N + n) x {N + n) matrices are inverse 
of each other 



-do)' \ -R-^DS-^ -R-^ 

Thus the system (4.4) is equivalent to the following one in normal form 

II J V -R-^DS-^ -R-^ Ms 



(4.5) 



(4.6) 



Since the second equation in (4.4) is obtained from the constraint equation by deriva- 
tion in t, we have that h{t, q, q) = const along the solutions of the first equation in (4.6). 
So, 6(t, q,q) =0 is satisfied provided it holds at the initial time to. Finally, we can get rid 
of the multipliers arriving at the following Cauchy problem which generalizes (2.10) 

q = F{t, g, q) 

(4.7) 

{q{to) , q{to)) = (go, go) e {(«, v) : b{to, u, v) = 0} 
where the first equation is the first equation in (4.6), briefly written. 

Example: the nonholonomic skate with rotational friction. We modify the 
Lagrangian function (2.11) so it has rotational dissipation 

L{t, X, y, 0, X, y,4') = \ (f + + 4>^) + x (4.8) 

where > is a new parameter. 

The nonholonomic equations are 

X — 1 — n sin 

y^-licos<t> , . 

+ A;0 = O ^ ' 

X sin (j) — y cos = 

Only the third equation is different from the one in (2.13). Instead of (2.14) we then have 

X — cos^ sin {x cos + y sin 0) 

y — cos sin + cos [x cos + y sin 0) -^^-^ 
= —k(j) 

(x(0), y(0), 0(0), ±(0), y(0), 0(0)) = (0, 0, 0, 0, 0, w) 
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The solution is now 



(4.11) 



m = f (1 - e-^') , x{t) = \ cos (0(r)) dr^ , 

y{t) = f sin (0(0) t cos(0(t)) drd^ . 
Jo Jo 

Notice from x{t) that the skate shdes eventually down. Indeed, if 

W f 7T "1 

lim 0(t) = — ^ i- + rmr :m eZ} then lim x(t) = +oo ; (4.12) 

t— >-+oo I 2 J t— >-+oo 

and otherwise x{t) converges to a strictly positive finite limit as t ^ +oo. 
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